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ABSTRACT 


For the case of homogeneous, isotropic magnetic field fluctuations, 
it is shown that most theories which are based on the quasi-linear and 
adiabatic approximations yield the same integral for the Fokker-Planck 
coefficient for the pitch angle scattering of cosmic rays. For example, 
despite apparent differences, the theories due to Jokipii and to Klimas 
and Sandri yield the same integral. It is also shown, however, that 
this integral in most cases has been evaluated incorrectly in the past. 
For small pitch angles, the errors in previous evaluations are for- 
tuitously of minor importance. For large pitch angles, however, these 
errors become more significant, and for pitch angles of 90 , the ac- 
tual Fokker-Planck coefficient contains a delta- function which has 
been overlooked in the past. The implications for these corrections 
on our ability to relate cosmic-ray diffusion coefficients to ob- 
served properties of the interplanetary magnetic field are discussed. 
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Introduction 


In recent years a number of seemingly different theories have 
been developed for relating the cosmic-ray diffusion tensor to ob- 
served properties of the interplanetary magnetic field (e,g* Jokipii, 
1966, 1967, 1968, 1971, 1972; Hall and Sturrock, 1967; Hasselmann and 
Wibberenz, 1968; Klimas and Sandri, 1971, 1973a, b). These theories 
are derived by using versions of both the quasi-linear and the adia- 
batic approximations o In the quasi-linear approximation the inter** 
action between the particles and the random magnetic field is assumed 
to be sufficiently weak so that to lowest order of interest a particle 
executes an unperturbed helix about the mean field while interacting 
with the random field. In the adiabatic approximation the cosmic-ray 
distribution function is assumed to vary on a time scale which is 
longer than and distinct from the time scale over which a particle 
experiences a significant interaction with the random magnetic field 
(see Appendix A to this paper for details concerning how these approxi- 
mations are used, and for conditions under which they are valid). 

In this paper we show that despite apparent differences there are 
substantial areas of agreement between the various theories. For ex- 
ample, we show that in the case of magnetic field fluctuations that are 
homogeneous and isotropic the theories due to Jokipii (1966) and to 
Klimas and Sandri (1971, 1973a, b) yield the same basic integral for 
calculating the Fokker-Planck coefficient for pitch angle scattering 
(< (Ap,)^ >/At , where is the cosine of a particle’s pitch angle). 

We also point out, however, that this integral for <(Ap)^>/At in 
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most cases has been evaluated incorrectly in the past. For example, 
there is no foimal justification for the approximations due to Jokipii 
(1966, 1971), which permit this integral to be expressed simply in 
terms of a one -dimensional power spectrum of magnetic field fluctua- 
tions, evaluated at a resonant wavenumber. By using both analytic and 
numerical techniques, we correctly evaluate <(A^,)®>/it for a general 
form of the power spectrum, which should provide a good fit to observed 
power spectra. In the limit of low rigidity, where the ratio of the 
correlation length for field fluctuations to the particle's gyro-radius 
is large,; the correct result for <(Ap,)^>/At differs from that due to 
Jokipii (1966) by a factor of [p,| (for p,^0). At all but very small 
values of this correction will introduce only a small numerical 
change in predicted values of <(iy,)®>/At. However, in certain schemes 
for averaging <(A|j,)2>/At over p,, and for forming the diffusion coefficient 

f 

for propagation parallel to the mean field (Jokipii, 1966, 1971; Earl, 
1973), this extra factor of ||i( introduces strong divergences as 
These divergences no doubt result from the inadequacies near ^ = 0 of 
the quasi-linear and adiabatic approximations. We also show that 
quasi-linear/adiabatic theory does not in general predict that 
< >/ At - 0 at |j, = 0 (e«g. Jokipii, 1966, 1971) or that <(Ajj,)^>/At 

diverges continuously to infinity as (Hasselmann and Wibberenz, 

1968). Rather, this theory predicts for physically realistic power 
spectra that < (A^,)^ >/At contains a delta- function, 6(|j,). Finally, 
we discuss how the above corrections affect our ability to determine 
an accurate diffusion coefficient for propagation parallel to the mean 
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field, and how they affect our ability to eliminate the inadequacies 
in current quasi-linear/adiahatic theory. 

Agreement Among the Theories 

The formula for < >/At due to Jokipii (1966) contains a 

minor algebraic mistake. In Appendix A to this paper we rederive the 
correct expression for < (A|j,)VAt that results from applying the quasi- 
linear and adiabatic approximations. We follow here procedures very 
similar to those used by Jokipii (1966). The resulting <(Ay,)^>/At 
for homogeneous, isotropic field fluctuations, which is the case con- 
sidered commonly by Jokipii (1966) and Klimas and Sandri (1971, 1973a, b), 
is given by 

< ^ J dX[cos(e\)[2b(r/Xc) +P^a(r/Xc)] - P®a(r/Xc)] (D 

At \c ° 

Here e = X^/tg> the ratio of the correlation length of the field fluc- 
tuations to the particle's gyro-radius; v is particle speed; and 

tP = B ^ /B ^ is the ratio of the mean square field strength to the 
T “rms' o 

mean field strength squared. The functions a and b of the scaler 
quantity r/X^ are the arbitrary functions that appear in the two-point 
correlation tensor for homogeneous, isotropic fluctuations (Batchelor, 
1956): 

Ci.(r) s <Bi(R)B^(R+r)>/7fB 3 

ij - - J (2) 

= a(r/Xc)rirj + b(r/Xc)6ij 
The functions a and b in (1) are to be evaluated at 
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^ j2 

r/Xc = + z (1 - cos(eX))] = [ji^X® + p^] (3) 

which is the linear distance from a point in a particle’s helical orbit 
along the mean field (chosen to be the origin) to the particle’s posi- 
tion at a dimensionless time X = Vt/X^ later. 

In the expression for < >/At due to Jokipii (cf . equation (18) 

in Jokipii (1966) or equation (32) in Jokipii (1971)) the sign of the 
second term containing a(r/X^) in (1) is incorrectly positive. Jokipii 
(1966, 1971) subsequently neglects these terms containing a(r/X^) so 
that the formulas for <(A[ 4 »)®>/At in terms of the one-dimensional 
power spectrum, which are normally attributed to Jokipii, are un- 
altered by the correction. However, as we discuss below, it is 
incorrect to neglect these terms containing a(r/X<.) at low rigidities, 
where e ^ 1. 

Klimas and Sandri (1971, 1973a, b) have also attempted to relate 
the diffusion tensor to observed field properties. Unlike Jokipii 
(1966, 1967, 1968, 1971, 1972), Klimas and Sandri (1971, 1973a, b) average 
first over pitch angle (on assuming that the particle distribution 
function can be expanded in terms of Legendre polynomials), and then 
they integrate with respect to time, i.e. they integrate over \ from 
0 to Qo. In Appendix B to this paper, we shdw by reversing the order 
of integration that Klimas and Sandri find a diffusion coefficient 
for propagation parallel to the mean field: 



K = r dp, r dxCcos(eX)[2b(r/Xc) + - p^a 

II 9 L-'o Xc ° 

(4) 

The first integral in (4), over X> can be identified with <(Ap,)®>/At 
given in (1). The second integral,, over p, 

(5) 

is identical with equation (4) in Jokipii (1968) or equation (46) in 
Jokipii (1971), and it performs an averaging of < (Ap)® >/At ' over pitch 
angle when the particle distribution function is expressed as an expan- 
sion in Legendre polynomials. The full expressions for K and 
< (Ap)^ >/At due to Klimas and Sandri (1971, 1973a, b) are then identical 
with the corresponding corrected results due to Jokipii (1966, 1968, 1971). 
However, Klimas and Sandri (1971, 1973a, b) have analyzed these integrals 
(e.g. (4)) in more detail than Jokipii (1966, 1971) and have been able 
to show formally that within the context of the quasi-linear approxi- 
mation a parallel diffusion coefficient exists and propagation parallel 
to the mean field can be approximated by a Markovian process. 
Approximations to the Integral for < (Ap)^ >/At 

The integral in (1) is quite complicated and therefore is probably 
not of much practical value unless it can be simplified by further 
approximations. At high energies, where € <<1 , simple approximations 
to (1) are possible and most theories agree that the resulting diffusion 
coefficient has a rigidity-squared dependence (e.g. Jokipii, 1966; 
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,1 , < (A;i)2 > 

I — 7 , — 

o At 



Klimas and Sandri, 1971). At low energies, where e>>l, valid 
approximations to (1) are more difficult to obtain. However, it is 
in this regime that it is most important to know <(A|j,)2>/At since 
in the interplanetary medium galactic cosmic rays that experience sig- 
nificant modulation and most solar cosmic rays should satisfy e>>l. 

Jokipii (1966) suggested that when e>>! the integrand in (1) 
could be expanded in powers of 1/e. To lowest order, then, the terms 
containing a(r/X^) were neglected and b(r/X^) was evaluated at r/X^, = [j,X 
or (1) becomes 


) 

The integral in (6) is simply proportional to the one -dimensional 
power spectrum of fluctuations in the component of the field perpen- 
dicular to the mean field direction, evaluated at the resonant wave 
number kX^ s/n* 


<(A|l)2> 2(1. 

||■■^|| I I ^ 

At 


■^i^)Tpe^v 




J dX'cos^^^^b(X* 


( 6 ) 


< (Aj.t)^>^ 2TT(l-4i^)7]^e^v 
At X^H,] 


Pj^(k=l/p,rg) 


(7) 


A similar expansion is implied in the work of Hall and Sturrock (1967) 
(cf. Equation 3.11 in their paper). 

This approximation technique, however, is unlikely to be valid 
except for certain special forms of the correlation function (not the 
isotropic case considered here). By neglecting terms containing 
a(r/Xc) in (1) and evaluating b(r/Xc) at r/X^ = )j,X, we are in effect 
making a guiding center approximation. We evaluate the correlation 
function only at the particle's guiding center, or equivalently we 
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assume that knowledge of the field at this point is sufficient to deter- 
mine the field as seen by the particle. This approximation should work, 
however, only if the field varies slowly on the scale of a particle’s 
gyro-radius. Yet in our final answer (6) and (7), <(Au,)^>/At is de- 
termined by a resonant interaction in which a particle responds only to 
fluctuations that have a scale-size on the order of a gyro-radius. 

We can see the difficulty with this approximation by examining 
(1). For large c, the cos(eX) in the term b (r/Xc)cos (eX) varies rapidly 
with \ between positive and negative values. Thus, as the integrand adds 
to form the integral, considerable cancellation will occur between these 
positive and negative parts. To determine <(A|^)^>/At, then, we must 
know the integrand to higher accuracy than results from performing an 
expansion in powers of 1/e and keeping only the first few terms. 
Mathematically this is a common problem where a valid asymptotic ex- 
pansion of the integrand does not necessarily lead to a valid 
asymptotic expansion of the integral, 

Jokipii (1972) has shown that there are special cases of the cor- 
relation function for which (6) and (7) are valid. If the field 
fluctuations vary effectively only with distance along the mean field 
direction, then the correlation length perpendicular to this direction 
is essentially infinite. In this case our objections to the guiding 
center approximation used in deriving (6) and (7) no longer hold, since 
now knowledge of the field at the guiding center is sufficient to 
specify the field as seen by the particle. It is debatable, however, 
as to whether actual interplanetary fluctuations have such an aniso- 
tropic correlation function. 
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In general, then, it does not appear that we can determine 
< >/At in the limit of large e by making a simple approximation to 

the integrand in (1), and we mast thus apply numerical and analytic 
techniques to a study of the complete integral* We use in these studies 
a fairly general form of the correlation function, which results in 
power spectra that are in rough agreement with observed spectra* We 
take 

1 / r \v-3/2 

%Xe^r(v4) (h;) ^-3/2(rAc) (8 

and 


.(^)K„.3/2(r/V)3 (9) 

where T(z) is the gamma function and K^(z) is a modified Bessel function 
of the third kind and index v* By applying the appropriate Fourier 
cosine transforms to (8) and (9) (see Batchelor, 1957), we find that 
the corresponding one -dimensional power spectrum for fluctuations in 
the field component parallel to the mean field direction is given by 


Pii(k) 


r(v)Xc 

3/^ r(v-i) 


(1 + (kXc)®)"'' 


( 10 ) 


and the corresponding one -dimensional power spectrum for fluctuations 
perpendicular to this direction is given by 


Pi(k) 


P||(k) 7 ^ 2v(kXc)^ \ 

2 \ 1 + (kXc)^/ 


( 11 ) 


We must choos.e v here to be greater than g so that there is finite 
energy in the field fluctuations (we could choose v and impose a 
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cutoff in the power spectra at high frequencies without altering 
our results). However, since we can choose v within the range, e.g, 
l/2<v<3/2 these oower spectra can be adjusted to give a good fit to 

observed spectra (cf. Jokipii and Coleman, 1968; Sari and Ness, 1969, 
1970). 

Using the a(r/X^) and bCr/Xj,) given in (8) and (9) respectively, 
we have integrated (1) numerically for a wide range of values of e, v 
and jj,. In the regime c > 1, we find that, to a good approximation (for 

li>0) 


< (Api)2> 2n(l-^j,3)'T|2e2v 


Pj_(k = 1/p.rg) 


where Pj^(k) is given in (11). That is, the actual <(A|j,)®>/At for this 
correlation function differs from <(Aj^)®>/At given in (7) (due to 
Jokipii (1966, 1971)) only by a factor of |p,j. In Figure 1 we have plotted 
the ratio of < (A^,)^ >/At that results from our numerical calculations to 
<(Ay,)®>/At that results from (12) versus for e in the range 
1 £ e ^ 10 and v in the range 3/4^vS5/4. As can be seen in this figure 
to a good approximation this ratio simply equals (x (g, 0). 

For this general correlation function, then, the answer for 
< (A|x)® >/At due to Jokipii (1966, 1971) is numerically not significantly 
different from the actual answer, at all except very small values of g,. 
It must be emphasized, however, that to a great extent the success of 


(7) is fortuitous. The function b(r/X^) evaluated at the full 
given in (3) is frequently an order of magnitude larger than b(r/X^) 
evaluated simply at the guiding center r/X^ = |iX. This result disproves 
the conjecture by Jokipii (1971) that higher order corrections to b(g^X) 
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can be neglected for v<l. The terms a(r/Xj,), which are neglected in 
Jokipii's approximation, are actually of the same order of magnitude as 
the full h(r/X^) terms, but are of opposite sign. The various terms in 
(1) simply add in a complicated manner to form the < (Ap,)® >/At which is 
given in (12). 

We have also studied (1) analytically. In Appendix C to this 
paper we derive the exact analytic expression for the < (Ap,)^ >/ At that 
results from using the correlation function given in (8) and (9). In 
the limit that this expression reduces to 


< (Ap,)5 > 

At 


1 + 3/2) ^ve^ ^ 1 

3 r(v - 1/2) Xc W ' 


where 



g(m,v,y) 


2ti= 


[g(m + l,v,y) + g(m - l,v,y)] 


(13) 


^ r(m + l/2)r(v + 1 - m) 

g(m,v,y) = y™ 

r(v + 3/2)r(2m + 1) 


+ 1/2; m - v, 2m + 1; y) 


(14) 

1^2 3/2; V + 2 + m, V + 2 - m; y) 

r(v+2-Hn) 

and y == tn? (1-^i^ ) . The function 3,Y; y) is a generalized hyper- 

geometric function* It can be seen in Appendix C that the sum over m in 
(13) is a sum over harmonics of the resonant wavenumber 

The e-dependence of (13) is particularly simple and has the same 
functional form as the e-dependence given in (12) (for large e). 
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Although we can verify (12) numerically only for l^e^lO, we can conclude 
by using (13) that (12) should hold for arbitrarily large values of e. 

The dependence of (13) on |j, and v is contained in the complicated sum- 
mation over hypergeometric functions; a summation that does not appear 
to simplify in any meaningful way. To the best of our knowledge, these 
generalized hypergeometric functions are not special cases which reduce 
to simpler functions, nor does it appear likely that the hypergeometric 
functions can be combined in any simple manner. For example the two 
hypergeometric functions which appear in g(m,v,y)(14) linear inde- 

pendent solutions to the same differential equation. Of course, we can 
always calculate the hypergeometric functions in (13) and (14) by using 
numerical techniques. Upon summing typically over the first ten har- 
monics, we can verify that in the range 3/4^v^5/4 and p, > 0, < (App)^ >/4t 
given in (13) is well approximated by < (Ap,)^ >/At given in (12). 

The extra factor of ip | in (12) as compared to (7) can have a sig- 
nificant effect on the pitch angle averaged diffusion coefficient in 
certain averaging schemes. Recently, Earl (1973) has argued that 
< (Ap)^ >/At should be averaged over pitch angle to form by use of 
the formula 



p^i’_(Hf)d^ 
o < >ll^t 


(15) 


rather than by use of (5). The averaging scheme in (15) was originally 
suggested by Jokipii (1966). However, as can be seen by substituting 
(12) into (15), for any v^l/2, the first integral in (15) will tend to 
diverge at its lower limit, p = 0. That is, IC|| as determined from (15) 
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does not appear to exist for any power spectra which have the general 
form given in (10) and (11) and which have finite energy in the field 
fluctuations. We should perhaps be skeptical then about this averaging 
technique. More probably, we should be skeptical about the ability of 
our theory to adequately treat the behavior of particles with pitch 
angles near 90® (p> 0). In quasi-linear theory, particles resonate 

only with fluctuations at wavenumber k = l/|j,r . As ^ ^ 0, the reso- 

o 

nant wavenumber tends to infinity, Pj^(k = 1/urg) -> 0 and thus 
<(Ajj,)^>/At ^ 0 (exactly at ix=0, <(Ay,)®>/At 0, as we shall discuss). 

However, if we were to improve quasi-linear theory with, for example, 
some resonance broadening technique (Duprey, 1966; Weinstock 1969; 

Jones et al., 1973; Volk, 1973) then we should expect that .< (Ap,)^ >/ At 
would remain non-zero as p.^0 and the divergence in (15) disappears. 

The Behavior of < (Au,)^ >/At at u =0 

Finally, we comment that theories which use the quasi-linear and 
adiabatic approximations for isotropic turbulence predict neither 
< (Ap,)^ >/At = 0 at |j, = 0(e.g. Jokipii, 1971), nor a continuous divergence 
in < (Ap)^ >/At for small p (Hasselmann and Wibberenz, 1968). In Ap- 
pendix D to this paper, we show that for an arbitrary isotropic 
correlation function, to lowest order 1/e, (1) predicts that 


< (Ap)^ > 

At 


TTlp V 
2 


no 

6(p) J dkkP|| (k) 
0 


at |j, - 0 


(16) 


where 6(p>) is a delta function. 

We can obtain an expression valid for all € by following the tech- 
niques outlined by Hasselmann and Wibberenz (1968) where the perturbation 
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magnetic field rather than the correlation tensor is expressed as 
a Fourier transform, Hasselmann and Wibberenz obtain an incorrect 
result because of an error when integrating over the wave -number 
phase angle (Volk, et al., 1973). In Appendix E to this paper we 
apply these techniques correctly, and obtain for |i = 0: 


<(AU,)®> TT(€T!)‘ 


At 




6(m,)v J dkP (k) 




(17) 


where 

^(x) sJi^Cx) - Jo(x) J 2 <x) 

and the ordinary Bessel functions of the first kind. 

We note that (16) and (17) do not describe a resonant interaction, but 
rather a particle responds here to fluctuations at all wavenumbers. As 
a particle spirals along its orbit it will of course experience gra- 
dients in the field due to all fluctuations, and it will have its 
pitch angle altered as a result of these gradients. In quasi-linear 
theory, however, the effects of these gradients at one point in a par- 
ticle’s orbit will in general be cancelled by the effects of the 
gradients at some other point further down the field, so that for |jt 0 
< ( A|j,)^ >/At 0 as At-^o®o Of course, the one exception to this cancel- 
lation comes from fluctuations at the resonant wavenumber, which can 
give rise to an accumulated change in the pitch angle, i.e. < (A|a,)^ >/At 
remains non-zero as At^^. For p, = 0, however, quasi-linear/adiabatic 
theory requires that a particle simply remains in one position along the 
mean field for infinite time. A particle, then, experiences the field 
gradients only at this one point, no cancellation of the type described 
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above will occur for any wavenumber fluctuation! and < >/ t,t tends 

to infinity according to (16) and (17). We note further that this delta- 
function at [ 4 , = 0 will not remove the divergences in (15) that occur in qua si- 

linear theory; rather this delta function makes the integral poorly 
defined. Also, this delta- function will not occur in Jokipii's special 
case (Jokipii, 1972) where the correlation length perpendicular to the 
mean field direction is infinite. A particle with p, = 0 experiences a net 
change in the random field component perpendicular to the mean field direction 
during its circular orbit, only if the random field along the mean field 
direction is changing with distance. That is, for a p = 0 particle to 
be affected by the random field on the average, must not be iden- 

tically equal to zero, which is the case in these infinite correlation 
length models. 

We do not wish to imply that this delta- function should necessarily 
be a part of a physically realistic derivation of < ( >/ At . However, 
the physics that gives rise to this delta function (i*e* changes in 
field magnitude, Pj|(kV0) must be taken into proper account. 

Concluding Remarks 

We have shown that despite apparent differences most theories which 
are based on the quasi-llnear and adiabatic approximations yield the same 
basic integral for <(A^,)^>/At for the case of homogeneous, isotropic 
field fluctuations. For example, the theories due to Jokipii (1966, 

1971) and to Klimas and Sandri (1971, 1973a, b) yield the same integral. 

We have also shown, however, that this integral in most cases has been 
evaluated incorrectly in the past. For values of near unity, these 
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previous evaluations fortuitously yield values of <(A^)®>/At close 
to the actual values. As however » the errors in previous evalua~ 

tions become more significant, and at |ji = 0 the actual < (A^)® >/At 
contains a delta-function which to date has been overlooked. 

We should of course be dubious as to whether the quasi- linear and 
adiabatic approximations can be applied for dealing with particles at 
linear zero. As a result, the corrected quasi-linear/adiabatic theory 
presented in this paper is not expected to be a physically realistic 
theory for all p, particles. With a small change in field strength small 
p, particles can reverse direction, in violation of the quasi-ltnear ap- 
proximation. These particles also move relatively slowly along a 
correlation length, and thus there is no clear separation between the 
time-scale for a 'scattering' of a particle and the relaxation time of 
the particle distribution function, as is required for the adiabatic 
approximation to hold. No doubt for particles with small p, quasi- 
linear/adiabatic theory should be modified, and replaced by a theory 
that treats the particle trajectories in a more realistic manner. 

Due caution, however, should be exercised in making any corrections 
to quasi-linear/adiabatic theory. We note in particular that none of 
the difficulties in evaluating <(Ap)^>/At that we have pointed out in 
this paper arise in models where the correlation length perpendicular to 
the mean field is essentially infinite (in at least one direction). Such 
models have been considered recently in some detail by, e.g. , Jokipii 
<1972, 1973), Jones, et al. (1973) , Vblk (1973) , and Volk, et al., (1973). 
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In any physically realistic model, however, where P||(k) is not 
identically zero, the corrections pointed out in this paper are expected 
to be important. In particular, quasi-linear/adiabatic theory will pre- 
dict a delta- function in < (Ap,)® >/At at p, = 0. Any modifications to 
quasi-linear/adiabatic theory that are developed by using one of these 
infinite correlation length models will be hard pressed to deal effec- 
tively with this delta-function. As a result, such modifications are of 
doubtful importance for the development of a physically realistic theory 
far the behavior of particles with n, near zero. 

We note also that current quasi-linear/adiabatic theory is not ex- 
pected to yield accurate pitch angle averaged diffusion coefficients 
(K^|) even though the inadequacies in the current theories may occur 
only in a narrow range of pitch angles around p, = 0. We saw above that 
the extra factor of p, that we found in < (Ap,)^ >/^t caused one of the 
averaging schemes for forming (cf. Equation (15)) to diverge as p-»0. 
It can also be shown that for reasonable power spectra, the other avail- 
able averaging scheme (Equation (5)) is dominated by the contribution 
from the delta-function, at low rigidities (Kliinas and Sandri, 1973c), 

In conclusion, then, no complete theory appears to be available 
at this time for accurately determining cosmic-ray diffusion coefficients 
from observed properties of the interplanetary magnetic field. 
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Appendix A 

Follcjwing Jokipii (1966), we calculate < >/At by using the 

Lorentz force equation 

^ ± (V, X B J (A-1) 

dt me 

Here v and v , are particle velocity parallel and perpendicular to the 

“'ll 

mean field direction, respectively; q is particle charge; c, the speed 
of light, and m, the relativistic particle mass. The magnetic field per- 
pendicular to the mean field is considered to have a random component 
which is denoted by and which has components and By in a Cartesian 
coordinate system with z-axis along the mean field. The average value 
of <B^>, is zero. (Angular brackets hereafter denote ensemble 
average . ) 

To make ready use of (A.l) we must make the so-called quasi-linear 
approximation. We assume that the interaction between the particles and 
the random field is sufficiently weak so that to the lowest order of 
interest a particle executes an unperturbed helix about the mean field 
while interacting with the random field. Then, v^ on the right side of 
(A.l) is given simply by 

Vj^ = Vj^(ey cosoflt - exSiniot) (A. 2) 

where e^ and ey denote unit vectors in the x and y directions, respec- 
tively; Vj^ is constant, independent of time; and, u) = qB^j/mc is the 
gyro- frequency of a particle in the mean field Bq* According to Klimas 
and Sandri (1973a) the quasi-linear approximation should hold for e^l 
provided that e^T|®<<l, where e is the ratio of the correlation length 
for field fluctuations (X^.) to a particle’s gyro-radius, and T\ is the 
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ratio of th© mean square field strength to the mean field strength 
squared^ At low rigidities, where 8^1, Klimas and Sandri (1973a) argue 
that the criterion for the validity of the quasi-linear approximation is 
only that if <<1, provided that the particles are of sufficiently small 
pitch angle so that they do not mirror. 

By using (A. 2), the mean square change in a particle's pitch angle 
which occurs in a time interval 0 to At can be obtained simply by inte- 
grating (A.l) twice, or 


< (Ap,)2 > 


<(Avi|)^> 

v2 


2-n2 it ^t-t 

= — — i r dt r dt'[cos lut cos (Dt'C (t,t’) 

+ (cos out sin U)t' + sin out cos out') Cjjy(t,t') 

+ sin out sin out' Cyy(t,t')] (A. 3) 

Here, we have replaced <3^ (t)Bj (t ' ) > with the corresponding compo- 
nents of the two-point correlation tensor for field fluctuations. 


(r(t,t')) 


<Bi(R)Bj(R + r(t,t'))> 

Bo® If 


(A.4) 


The field fluctuations are assumed to be homogeneous, so that is 

independent of some arbitrary pbsition vector R and a function only 
of r(t,t'), the vector separating a particle's position at time t and 
at t ' . 

In the next step in our derivation we change variables in (A. 3) 
from t and t ' to t and t = t'-t, and then we reverse the order of inte- 
^fter a considerable amount of tedious but straightforward 
algebra, we find that (A. 3) becomes 



V./ 


< = 


pLt 

— .L 


r />At-T r 

a ‘‘4' 


cos CUT 


(^xx ■*■ *^yy) 


+ ( 
+ ( 


cos CUT COS 2cot - sin cur sin 2cut) 




XX 


sin (UT cos 2(ut + cos (Ut sin 2W)t) C 


■b 

xyj 



r 


cos (JUT 


(Cxx + ^yy) 


1 / /s . • V XX ^yy 

+ (cos CUT cos 2cut + sin cur sin 2cut)| — — 


+ (cos CUT sin 2tut - sin cut cos 2cut) C 


l ^xx - Sy j 

■]} 


If we further assume that the fluctuations are isotropic, then 
Cij (r) = b(r) 6ij + a(r)rirj 

where a(r) and b(r) are arbitrary even functions of the scaler 
separation r, which in turn is a function only of t given by 


(A.5) 


(A. 6) 


2v 

r3(T) = V ^ (1 - cos cut) 

11 0)3 

= V|®t 2 + p3(T)X^2 (A. 7) 

The first integral in (A.5), then, can be done exactly; it is inde- 
pendent of a(r(T)) and b(r(T)) and involves integrating only simple 
trigonometric functions. The resulting terms are either linear in 
(At - t) or proportional to 1/cu. These latter terms can be neglected and 
we find that 


< (Am,)2 > 


cu^v.^tP 

— 


(At - t)[cos cuT[a(r(T))P^X^.+2b(r(T))] - a (r (t) ) ] (A. 8) 
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Finally, we make the so-cal led adiabatic approximation* The 
integrand in (A. 8) is assumed to fall-off with a characteristic time 
T, so that for At>>r, we can approximate At-T as At and extend the 
upper limit of integration to infinity. In terms of a dimensionless 
time X = vt/Xq^} (A. 8) then reduces to (1). 
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Appendix B 

In this appendix we show that the result of Klixnas and Sandri 
(1971, 1973a, b) for K|| can be expressed in the form given in equation 
(4) of the text. 

Upon replacing Klimas and Sandri 's characteristic length \p with 
our Xp (\p = /rr X^,/6), we see from equation B.4 in Klimas and Sandri 
(1973b) that these authors find that 


iC = i ^ ^ ^ 

II 3 Tf 


(B.l) 


By using equation B.3 in this same paper, and equation 3.19 in Klimas 
and Sandri (1971), we find that 




A = JdX 2Jj_(e,X) cos (eX) - 8X2Ii(€,X) 


'cos(eX) - 11 
~eX 


(B.2) 


where, by using equations (3.22), (3.24), (3.15), and (3.16) in this 
latter paper. 


3 

~ J^d|J. fa(r/X(,) (l-y«^) 


(B.3) 


and 


Ii = I J dij, a(r/X )(l-u®)' 
o o ^ 


(B.4) 


We note here that in Klimas and Sandri 's notation, our a(r/X(.) and 
b(r/Xj.) ace B(r/X^) and A(r/X^.), respectively. Finally by substituting 
(B.2), (B.3) and (B.4), from above, into (B.l), and upon reversing the 
order of integration, we obtain (4). 
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Appendix C 

In this appendix we evaluate analytically <(AiJ.)2>/At when the 
povjer spectrum of magnetic field irregularities is given by eqs. (10) 
and (11). Starting from (1) we have 


< (A^,)3 >/ At 


(l-|i®)Tf e^v 


(C.l) 


where 


I = r dX {coseX[2b(r/Xc) + a(r/Xc)p^l- a(r/Xc)p®} 
*^0 


and 


r 







1 

+ p^)! 


(C,2)^ 


(C.3) 


Because r, a and b are such complicated functions of (C*2) cannot 
be integrated direct Ly* One possible approach would be to transform 
a and b into three-dimensional Fourier- transforms of the power spectra. 
The functional dependence of I on r would be contained in terms propor- 
tional to exp i k*r and no terms would be functions of jx[* However, 
because it is not possible to write the last term in (C,2) in terms of 
the correlation tensor, we proceed by first transforming I into a new 
coordinate system in which the magnitude of remains unchanged, but 


the components are 
rz ^y 

■-• = nX, — = p, 7- = 0 

Kc Ac 


(C.4) 


We can now define an isotropic correlation tensor j (x) with com- 
ponent s 

Rij(r) =b(r/Xc)6ij + a(r/X^)rirj 


(C.5) 
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so that 

Rxx = b(r/\^) + p^a(r/Xp) (C.6) 

Ryy = b(r/Ac> 

If we identify the scalar functions a and b in (C.5-6) with a and 
b in (C.2), we have 

I = ;“dX[cos(eX)(R^^ + Ryy) - - Ryy)] (C.7) 

The components of the tensor Rij(r) can be expressed as three-dimensional 
Fourier- transforms of a power spectral tensor P. j : 

= KfiJ® e -- (C.8) 

SO that 


I = i / dX {e CpkCPxx + Pyy)e^-’-] - Jdk(P^j^ - p. 


ik*r- 


yy 


)e - -} 


(C.9) 


We now interchange the order of inte^gration, and write 

2|j, ' 

k‘r = k^XjjlxX + — k^X^ sin(eX/2) (C.IO) 

where n’ = (l-|i^)^- The integration over X can now be performed by making 
use of the following identity 


eiz sin <P = S ei“Pjm(z) 

m=-oo 


where Jjj,(z) are the Bessel function of the first kind. 


Therefore (C.9) becomes 

^ M + fyy’ - ea-m/2)] 


(C.ll) 


(C-12) 
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Note that and Pyy are even functions of so that the sum in 
(C.12) is non-zero only for even m. Therefore, after integrating over 
k^, the integral can be written as 

^ +00 +OT (C « 13) 

^ a. My ^ 

m=-“ -00 -® = — 


In general, the Pxx and Pyy are given by (Batchelor, 1956) 


P. . (k) = (k^6 - . - k.k .) 

ij 4nk^ iJ 1 J 


(C.14) 


where 


E(k) = k^ 


[dk \k dk / J 


(C.15) 


and P|| is given by eq. (10). One should note that this definition of P|| 

ensures that a(r/\ ) and b(r/X ) are given by eqs, (8 ) and (9 ). The inte- 
c c 

gral over ky can now be performed with the result 

4r(v+3/2)e2\ 

- r m3 I 

o 


I = 


3^(v-l/2)^L^l' 


^ Sm3 J dk^ (1+ + 


m-o 


}Jb ^ 6 


(ca6) 


where we have rearranged the summation so that it extends only over 
m^O. For non-zero m*, the m=0 term is zero* That this term is in fact 
a delta- function is shown in Appendix D. (For non-zero m, each term in 
the sum equals zero for ^ = 0.) We now restrict ourselves to and 

drop the m=0 term* The remaining integral over can also be performed 
(Gradshteyn and Ryzhik, 1965) so that 
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< (i|i)2 >/it 


4r(v+3/2)11^ve^ 
3r(v-l/2) \c 



(1 + 


2 \-v- 




g(m,v,y) 


+ 


(C.17) 

where g(m,Vjy) is defined by eq. (14) with 

y = (1 + m2e2/|^2). (C.18) 

In the limit of low rigidity (e>>l), (C.17) reduces to (13) in 
the text. 


2u® 


[g(m + l,v,y) + g(m - l,v,y)] 
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Appendix D 

We can derive (16) by expanding — — in (1), in powers of 1/e; 

i.e,, in powers of tlie gyroradius. We have demonstrated that this pro- 
cedure does not produce a rapidly converging series for the resonant 
interaction, but we will see here that a true power series in (1/e)® 
can be obtained for the delta -function term. As a demonstration, the 
first correction to (16) will also be constructed. In Appendix E we 
derive this delta function term from a different point of view for all 


e . 

The two functions, and b(“ ), which appear in (1) for 

Xc A.C 

^ are not independent of each other. From V = it 

At 

follows that 


CO r 

^ = 2 r dxx a (x) - (— ) a(~) 

Xc Wx. Xc 


(D-1) 


which we substitute into (1) to give 


11^ . , ^2 P.. ,(x) - (i) a(^)^ . 


-i (cos eX - 1)® (1 - p,®) a(— )] 


(D-2) 


We now construct the Taylor Series expansions of a(— ) and 
I(— ) = r dxx a(x) 

V r/X. 


(D-3) 


in powers of r^^ = "T ‘ results of these expansions 


are , 



no 


I (t^) = IW - i r a(^iX) + 

Ac 

and. 


a(-L) = a(M-X) +i r /(iltti) a'(^i\) + 


)j,X 

+ i -^o* ta”(nu 


a ' (m.X) 


] + 


M,X ' ' '" ' M.X 

These expansions are uniform in p, and X for 0 ^ 
Upon substituting these expansions into (D-2) we 
the series, 

<(Ay)2> _ <(A,i)S>[ -^1^2 <(A^)3>| 

, , • 1 + I + K—) 

At At e At Ig e 

It should be noted that the coefficients in this 
of e; as it stands at this point, (D-6) is not a 
expansion in 1/e. 

The leading term in (D-6) is given by, 


<(Ali)2>! „(l-ii^)TPe2v ” ,eX, 

A. - 2 J dX cos ( — ) b(X) , 

'o K n W 


or, in terms ofPj^(k), 


< (A|^)s> 
At 


(l-p,^)Tps^v 1 

= TT p (k = ) 

tirg 


M- 


(D-4) 


si, and 0 < X^“. 
are able to construct 

^ > I + . .. (D-6) 

At 

expansion are functions 
simple power series 


(D-7) 


(D-8) 


These are alternate expressions for Jokipii's (1966, 1971) approxima- 
tion to the resonant interaction given in (6) and (7). 
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The next term in (D-6) is given by. 


< > 
At 



T 

(1 y, )T1^e V j> [(a(n,X) + nX a'(n.X) ) + (D-9) 

T-*<» ° 

2 cos(eXX2 a(y,X) + y,X a'(|j,X) )+cos(2eXX3 a (]j,X) + u,X a ’ (M-X) )] 


The limit on T is included here to ensure that the change of integration 
variable, from X to r = n,X, which we are about to introduce, is carried 
out correctly. The last two terms in the integrand of (D-9) contribute 
to the resonant interaction discussed fully in the main text. Thus we 
ignore these terms here, and study the first term which represents the 
leading contribution, in 1/e, to the delta function term. This contri- 
bution is given by. 


< (Am,)® > 


At 


(l-^i2)Tl2e3v UT d 

• lim dr — (r a(r) ) 
dr 


M- X. 


(D-10) 


or 


<(A^)2> (l-^^2)Tfc®v 


At 


lim T a(|j,T) 
T 


(D-ll) 


If a(r)-»0 as r-»>», faster than 1/r, then, for |i = 0, X:m^T a(n,T) ^a(O) 

lim T = -. But, for ^^0, lim T a(^iT) = 0. Furthermore, 

->00 T 


^ ^ CO 

r du lim T a(^iT) = lim f ^(r) = J dr a(r) 

T^®o o 


(D-12) 


Thus, as a function of p,, lim T a(MlT) is zero for all w>0, is infinite 
at |j, = 0, and can be integrated to give a finite, non-zero number, i.e. 
it exhibits all the properties of a delta function, or 
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Tfe^v 

Ac U - ‘(l^> { ^(r) 

In terms of the power spectrum, 

< (A(i)^ > I nTpe^v 

— — |r — — «(n)/dk kp (k) 

o 


(D-13) 


(D-14) 


If we once again drop the contributions to the resonant interaction 
in the fourth order term in (D-6), we find. 


<(4u,)3> 

At 


(l-u^)^TPe2v , . 

‘ / Is, X «a"W) 


(D-15) 


For homogeneous turbulence in continuous fields, we must have a’(0) = 0 
(Batchelor, 1956), Therefore, 


<(Ap.)3> (i.,^3)3'f|3e3v a'(iiT) 

At , ■ T-« • »-i6) 

Finally, using reasoning similar to that presented above for the 
leading contribution to the delta function term, we find, 


< > I 

At I 4 


^e^v a’(r) 

6(p,) J dr 

Xc or 


(D-17) 


Since a'(r) is proportional to r for small r, the integral in (D-17) is 
finite. Equation (D-17) gives the X)((— ) ) term in a power series expan- 
sion of (D-6), Rather than continue this calculation to higher orders 
we present in Appendix E a different procedure which gives the delta 
function term for all e* 
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Appendix E 

We calculate below an alternative expression for the Fokker-Planck 
pitch angle scattering coefficient for isotropic turbulence and show 
that at p, = 0 it is proportional to a delta -function to all orders in e. 
We again begin with the equation of motion (A*l) and Fourier transform 
the random component of the magnetic field. One again uses the iden- 
tity (C.ll) so that 


< > 


4tt 


T T / 

f dX / dX' 


o o 


E(k) 

k2 


00 

E 

n=-“> 


00 ik||X£.li.(X-X ) 
E e 


iXcM-' ^ j ^kj^XcU.' 

[cose(X-X') - sin(\lf-eX)sin(il(-eX’)] (E.l) 

where E(k) is given by (C.15), and \|( is the wave-number 

phase angle. The integrals over ijr, X, and X' can be readily performed 
using the standard adiabatic approximation (cf. Appendix A) with the 
result 


jie(nX-raX') g-i\|r(n-m) j 


-e 


<(An)^> 

At 



fdk^k^ /”dk„ 


E(k) 

kS 


+ 00 


2 6(nk|,Xc+ me) 

m=-oo It 



(E.2) 


For n. 0 one can readily integrate over k^ . However, we want to in- 
vestigate (E.2) when y, = 0„ In this case we transform to polar 
coordinates and note that only the m=0 term in the summation is non- 
zero at jj, = 0. One can integrate over polar angle, leaving 
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< (A|i,)^ > _ ne^Tf V 
At 


2Xc^ 


E(k) ^/kXc\ 
6(^i) r dk-^ J^^[—) 
o k L \ e / 




Using the definition of E(k) one arrives at (17). 



Figure Caption - 


The ratio of the exact numerical integration of (1) ( < (fiji,)® > /At) 
to the approximate result (7) due to Jokipii ( < (A| 4 ,)® > /At) j is plotted 
against 0) for various values of e and v* To a good approximation 
the exact result is smaller than the approximate result by a factor of 









